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CURRENT DISTRIBUTION IN AN ARRAY OF SUPERCONDUCTING THIN FILMS* 


by I. D. Skurnick, J. H. Simmons, and A. R. Sass^ 

Lewis Research Center 

SUMMARY 

The conservation- of- fluxoid approach is used to derive the integral re- 
lations describing the current distribution in a general array of cylindrical 
London superconductors. This approach highlights properties of the kernel of 
the integral equation that permit the development of analytic solutions for a 
2 by 2 array of closely spaced thin films of rectangular cross section. It was 
found that the current is distributed across the width of the film in a non- 
linear , and nonsymmetric fashion, peaking at both edges. This peaking, in turn, 
dictates the switching characteristics of the system. It is shown that the de- 
pendence of the peaking on the separation distance between coplanar films can 
cause considerable variations in the critical current of the system I c as this 
separation is varied from 1 to 20 times the film thickness. 


INTRODUCTION 

A major problem in the design of thin film superconductive devices has been 
the difficulty in predicting the maximum current that can be carried by the 
system before it switches to the normal state. Studies of the critical current 
in single superconducting thin films revealed that the problem is complicated by 
a nonlinear peaking of the current density at the film edges (refs. 1, 2, and 3). 
Recently Sass and Skurnick, employing the concept of fluxoid conservation, de- 
rived the inhomogeneous Fredholm integral equation of the second kind, which 
describes the current density in both single and double cylindrical supercon- 
ductors (ref. 4) . The properties of the Fredholm kernel, in turn, permitted 
them to obtain a closed-form solution to the double-film case, a superconducting 
thin-film strip transmission line. In the strip line, as in the single film, 
the current density varied in a distinctly nonlinear fashion near the edges of 
the system. With respect to its widthwise variation, the current distribution 
in both cases is symmetric with respect to the center of each film. In this 
report the fluxoid conservation approach is used to extend the integral relations 
for the current distribution to a general array of cylindrical superconductors. 


A portion of this report was presented at the Third International Confer- 
ence on Magnetics, Washington, D. C. , April 21-23, 1965. 

*B. C.A. Laboratories, Princeton, New Jersey. 



The behavior of the Eredholm kernel permits the derivation of closed-form solu- 
tions for many cases of interest, in particular, the instructive example of a 
2 by 2 matrix of thin films. The solution indicates, signific ant ly, that the 
current distribution is strongly dependent on the proximity of neighboring films 
provided that such separations are on the order of 20 film thicknesses or less. 
This, in turn, suggests that the results obtained can be applied with reasonable 
accuracy to the more general case of a 2 by n matrix of thin films. 


SYMBOLS 


[The rationalized meter -kilogram-second system of units is used herein.] 
a unit vector (treated as dimensionless quantity) 

C constant determined by total current 

d thickness of film 

E electric field intensity 

l! magnetic field intensity 

I total current in individual film 

J current per unit area in x-y plane 

K(u,u*) kernel of one -dimensional Eredholm integral equation 


k Z/d 

l distance between centers of two films in strip-line pair 

m distance between centers of coplanar films 

P ln|(x - x') 2 Hb (y - y')2| 

q - 1 separation parameter, q = m/w 
r field point in x-y plane 

r 1 source point in x-y plane 

u dimensionless field point in one -dimensional analysis, (2/w )y 

u* dimensionless source point in one -dimensional analysis, (2/w)y T 

w width of film 

x,y Cartesian coordinates equivalent to r, when primed equivalent to r ! 
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p- 1 

London penetration depth 


k ( u ) 

coupling factor, J* 

K(u,u’ 

1 )du* 

A 

(Pd) 2 (w/d)(l/8it) 




permeability of free 

space , 

4jtXlO -7 H/m 

Subscripts : 



a 

array a 



b 

array b 



z 

z-component of vector 




GENERAL ARRAY 

-Consider the array of parallel cylindrical superconductors shown in fig- 
ure 1. In the following analysis the cross section of each superconductor is 



Figure 1. - General array of infinitely long cylindrical superconductors. Direc- 
tion of current flow in cylinders is arbitrary. (One cylinder has been enlarged 
to show typical contour of integration. ) 
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assumed arbitrary. If a direct current or low-frequency alternating current is 
p imp ed through the conductors, a quasi -static analysis is valid. The current is 
assumed to he in the z-direction; thus 

J = a z J z (x,y) 

Since the London equation ]f = q o P“ 2 (dj/dt) , which relates the electric 

field II, the current density J, and the London penetration depth 3 , is 

valid in each of the superconducting cylinders, application of Faraday’s induc- 
tion law to the dotted contour shown in figure 1 yields 


_d_ 

dt 


J(?) - j(o) - p 2 



(a z X dr") 


= 0 


( 1 ) 


where a z is the unit vector in the z-direction and dr" is a differential 
vector line element of integration. The quantity contained within the brackets 
in equation (l) is the London fluxoid; the vanishing of its time derivative ex- 
presses the principle of fluxoid conservation. Starting from zero field initial 
conditions requires the fluxoid associated with any contour contained within the 
superconductor to vanish. Thus, 

— > 

J(?) = J(0) + P 2 f H(r") • (a z X dr") (2) 

Jo 


The magnetic field at a point r" is readily determined by summing over the 
contributions made by the current elements in all the cylinders. Hence, 



Figure 2. - Four superconducting thin films with current directions indicated for 
arrays a and b. In array a the current enters system through films 1 and 2 and 
leaves through films 3 and 4; in array b the current enters through films 1 and 4 
and leaves through films 2 and 3. 
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( 3 ) 


H(?") =jj J(r') 
( all space ) 


G a x(?- -?•)] 

dr * 

2jt|?" - ?' | 2 


Substituting equation (3) into equation (2) leads to 


J(? ) = ° + £ |d 2? ' ( 4 ) 

(all space) 

where C is a constant determined by the total current. 

From equation (2) it is clear that the nonuniformity of the current density 
in the cylinders is due to the net magnetic flux per unit length through areas 
hounded by contours of the type described in figure 1. It should be expected 
then that if this flux is small; the current density would exhibit only mild 
variations in the cylinders. 


FOUR-FILM ARRAY 


The particular array shown in figure 2 depicts a 2 by 2 matrix of super- 
conducting thin films of rectangular cross section. The array can also be 
regarded as either two strip transmission lines in close proximity or two co- 
planar thin films deposited above a superconducting ground plane. Each film has 
a thickness d and width w. The spacing between the coplanar films is m - w, 
and the distance between the centers of a film and its image in a strip-line 
pair is l • The analysis is made for film parameters in a range of practical 
interest) that is, d/w < 0.01 and (Bd < 1. 


Applying equation (4) to the four film array yields 


J(x,y) = C + £- 
4Jt 


/ ' w/2 rd/2 
w/2 J- d/2 


J ± (x' ,y')F dxl dy 1 


rm+(w/2) /*d/2 

/ / J 2 (x* ,y' )P dx T dy 1 

J m-(w/2) d/2 


-(w/2) ./-d/2 


*/ 2 r(d/2 )-Z 

+ I / J 3 (x , ;y , )P dx 1 dy' 


r r 

y-w/2 y-(d/2)-Z 


m+(w/2) /*(d/2)-Z 


./a- (w/2) J-( 


/2) y-(d/2) - I 


J (x* ,y' )P dx' dy' 


( 5 ) 
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where P = ln|(x - x ' ) 2 + (y - y') 2 |. Figure 2 shows that the arrays a and 
h possess the following symmetry properties: 


Array a J(x,y) = -J(-x - Z,y) = j(x,m - y) = -<T(-x - l ,m - y) 


b 


J(x,y) = -J(-x - Z,y) = -J(x^m - y) = j(-x - l , m - y) 


( 6 ) 


Equation (5) can thus be greatly simplified by the following changes of vari- 
ables : 


Film 2 

3 

4 


m - y' 


"A 


x' -* -(x 1 + Z) 
both y* -> m - y' and x 1 -> -(x 1 + Z) 




(?) 


Incorporating expressions (6) and (7) into equation (5) leads to a more compact 
description of the current density in the films: 


2 r v/ 2 r d/2 
Ja,b(*,y> - C + - / / 


-w/2 -d/2 


X In. 


(x - x*) 2 + (y - y ') 5 


(_(x + x' + Z ) 2 + (y - y*) 2 J 


F. 


a,b 


v dx ' dy ' 


(8a) 


where 


F a = 


(x - x 1 ) 2 + (y + y' - m) 2 
(x + x* + Z) 2 + (y + y' - m)‘ 


and F, = — 
b F„ 


(8b) 


i ”1 


Since d < , it will be assumed for simplicity that J does not vary 

significantly with x. This assumption is not crucial and will be discussed in 
more detail at the end of this section. There is thus no loss of generality in 
setting x = 0. Integrating with respect to x', letting A = (pd) 2 (w/d) (l/8it) , 
k = Z/d, q = m/w, and changing variables such that u = 2y/w yield 


J a,b( u ) = C + A J J(u’ )K a ^’b( u ^ u * )du' 


(9a) 


where 
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» 

tf 


- . u. 

( 2k - 1) w 

u + u* - 2q 


, , ' ' W , n W 

+ tan - 1 - : = tan"- 1 - — — = tt 


(2k + 1) — 
w 

u + u* - 2q 


where the plus sign corresponds to array a and the minus sign to array b. 


Clearly, current elements in all four films contribute to the total mag- 
netic flux through any individual film. Figure 3 shows, however, that contri- 
butions from antisymmetric current elements in thin films tend to cancel them- 



Figure 3. - Contributions to magnetic field at P from distant antisymmetric source elements. (Film 
widths are not drawn to scale.) 
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selves out, particularly if these 
elements are distant from the 
field point under consideration. 
Let the currents at points S and 
S* he equal in magnitude but op- 
posite in direction. If Spg is 

the magnetic field at P due to 
the current at S, Hpgt the 

field at P due to the current at 
S * , and t^ > 7tg > it follows 
that 

jHp S |-|Hps'! 1Q 

I^PS 


(a) Absolute value of kerne! plotted for arrays a and b as function 
of u-u' for d/w = 0.001, k = 2.00, and q = 1.001. (Peak 
value has been normalized to unity.) 



u 


(b) Absolute value of coupling factor plotted against u for arrays a 
and b. (Parameters of system are the same as those used in 
fig. 4(a).) 

Figure 4. - Kernel and coupling factor. (Abscissa not drawn to scale. ) 


and 0.989 < cos 0 < l.OOO. Thus 
the great bulk of the flux through 
areas bounded by contours of the 
type shown in figure 1 is contri- 
buted by current elements compara- 
tively close to the field point. 

It should be expected then that 
the kernel would be sharply peaked 
about u = u 1 and decrease with 
increasing |u - u f | . Figure 4(a) 
shows representative one- 
dimensional kernels plotted as a 
function of u-u*. Although the 
kernels are not drawn to scale, 
the general features are evident 
and are supported by numerical 
analysis. 


Since l > d and m > w it can be shown that the one -dimensional kernels 
in equation (9b) are negative definite. Thus, equation (9a) can be rewritten as 


J(u) = C - A 



J(u* ) |K(u,u* ) |du' 


( 10 ) 


Equation (10) can be solved by the method of successive approximations. Starting 
with Jq(u) =0 it can easily be shown that 


J n (u) = C 


1 + 


n r l /*1 /*1 

( ’ A)m y 1 J-l * * * jfi l K(u m-l ,U m ) l * * • l K ( u ; u l)! du m * • * du l 


(Ha) 


* 


« 


The current density is given by 
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J(u) = lim J n (u) 
n-»°° 


(lib) 


provided that the series in equation (lla) converges as n -*• ». Thus, to the 
zeroth order in A the current density is assumed to be constant in the fi lm s. 
The first order correction considers the interactions of all the source points 
in the films with a particular field point. The second order correction con- 
siders the effect of all the source points on an undisturbed source point, be- 
fore the effect on a field point is computed. The higher order terms (in A) 
are further corrections for the interaction between current elements. 


It is helpful to evaluate the coupling factor defined by 

1 

'-1 

If*om equations (9b) and (12a) the coupling factor can be written as 


r 1 

J K(u,uq)du^ = k(u) 


(12a) 


K a,b( u ) = ( A 1 1 A l) 


where 


A x = In - 


Q=1 

+ 2(A 2 ± A£) 

(q + (u + l) 2 ] 


+ f k " ^V A i 1 A 7 




= 2 k-lf + ^ 1 


Q=2k+1 




u+1 


+ (a 2 ± a 2 ) 


Q=2k-1 


- ( A 2 ± A 2 ) 


(q if + (u - 1)2_ 


U^T - 4 + 2 l^ W )l tan_1 “T - tan_1 a 

w w 


Q=2k+1 


u - 1 


(12b) 


(13a) 


7 A 2 

“|u-2q+l 

L(« ;) + 

- 2q + l) 2 

r 7 2 

*iu-2q-l 

IQ i) + (u 

- 2q - 1)2 

\ -w/ 



- 4 + 



u - 2c[ + 1 



- tan 


-1 u - 2q - 1 



(13b) 
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A _ (1 + u 2) 1 2 _ i 

A 2 " ' i u J 2 d 2 w ' 


{[* * (tt) 


tan 


1 j; ^ 

«s 


i + 


1 - u 


Q - 

W 7 -I 


tan' 


Tq h a* } 


(13c) 


and 


«)« = -%. a. 


2 -, 


1 + 


2q - u + 1 


Q - 

V 


tan' 1 . 2 <L - . + 1 


Q - 

v 


1 + 


feq - u - 1 

v 


tan 


■l/2q. - u - 1) 2 


Q 


w 


i} 


- — !► (^) 

Q 


The plus sign refers to K a (u) and the minus sign refers to /^(u). 


Because of the behavior of the kernel, k(u) must also be negative definite. 
Since the kernel has a rather narrow effective width (fig. 4(a)), k( u) should 
remain approximately constant over the central portion of the films. On the 
other hand, when the distance to the edge of the film is on the order of the 
width of the kernel, k(u) is a strong function of u. This Implies that for all 
practical purposes field points in the central region of the film view the film 
as being infinitely wide. Figure 4(b) is a sketch of the absolute value of the 
coupling factor for typical values of the parameters of the system. Since 
k(u) csj k(0) over most of the width of the film, the following approximation sug- 
gests itself 

K ( Ul )| ^(ujupldu! =* |k(0) | |k(u) | (14) 



Thus in general 

•1 rl 


f±f-l ' f-1 ' K(Um ‘ 1 ' Um) '' * * l K ( U ' U l>l d %* * -|K(0)| m - 1 |K(u)| 


Combining equations (15) 


and (lib) yields 


J(u) 



A|k(u) 


” 00 

i L(-imb(°) 

_m=0 



(15) 


(16) 
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If A|k:( 0) j< 1 the infinite series in 
equation (16) converges absolutely, and 
the current density can he described in 
closed form by 


J(u) = C 


am u) 


1 + 7v|#c(0) | J 


(17) 


When ?\ | /c(0 ) | > 1, the series in equa- 
tion (16) does not converge, hence the 
method of successive approximations is not 
satisfactory. Figure 5 illustrates the 
range of convergence for equation (16) for 
d/w < 0.001. 


Figure 5. - Range of convergence for equation (16) for 
d/w <0.001. In this limit k c is independent of d/w. 
Permissible values of k are given by k<k c = 
^/(pd) Z ]+ 0. 25. 


A factor of particular interest and 
importance is the relative peaking of the 
current density at the edges of the film, 
that is, j(±l)/j(o). In contrast to the 
single film and strip line, it is clear 
that the distribution in arrays a and b can- 
not be symmetric. Considering the quali- 
tative features of the current distribution in the strip superconductor indicates 
that if the current in films 1 and 2 are in the same direction the greatest 
peaking should be expected on the outside edges of the system. On the other 
hand, if the current in films 1 and 2 run opposite to each other, the greatest 
peaking should occur on the inside edges. Figure 6 is a representative sketch 
of J(u) against u for both cases. 


From equation (17) the relative peaking at the edges is given by 

J(±1)/J(0) = 1 + a[|k(0)| - (k(±1)|] (18) 

In the limit where d/w < 0.001 it is possible to express the peaking in a form 
that explicitly demonstrates its dependence on various parameters of the system. 
From equations (12b), (13), and (18) the peaking on the outside edges can be 
shown to be 




1 + |a|k(0)| 



On the inside edges of the system equation (18) can he approximated by 


(19) 




« 1 + 


IMI 


i\ (e*) 2 * 2 1 r 1 1 

V 4/ 4« |_q(q -1) J 


( 20 ) 


where the minus sign is for array a and the plus sign is for array b. It is 
apparent that the dependence of the peaking on the spacing between coplanar 
films is only important when q - 1 is on the order of 20(d/w) or less. For 
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u 


Figure 6 0 - Normalized current density plotted against u for d/w = 
0.001, pd = 1, k = 2.00, and q = 1.005. (Edge regions of u-axis 
are not drawn to scale, > 


eter q - 1 in array ft than in array a. 
system the variation of J in the direc 


all practical purposes, the peak- 
ing at the outside edges of the 
system is unaffected fty the separa- 
tion between the coplanar films. 
Equation (20) is an excellent ap- 
proximation to equation (18) (i.e., 
it is in error fty less than 1 per- 
cent) for q - 1 > 5 (d/w). For 
smaller separations higher order 
terms in d/w must be taken into 
consideration. 

Since a superconductor starts 
to switch to the normal state when 
J exceeds a certain critical value, 
equations (19) and (20) and fig- 
ure 6 indicate that the critical 
total current will show a stronger 
dependence on the separation param- 
In calculating the total current in the 
.on of the film thickness should be con- 


sidered. If i(^) is the critical total current in array ft, calculated fty assum- 
ing a uniform widthwise current density distribution in the films (no peaking at 

the edges), it is obvious that iO 3 ) < I^j 3 ) . From equation (20) it follows that 


M 

c 

cu 


/l + (f 3 ^ 2 ! 

/ , V (pd) 2 k 2 i 

fk - 1 U w 

i 

1 4 

\ 4 / 4 at 

5(1 - 1 ) 


iT 


( 21 ) 


Figure 7, plotted fty combining equations (13) and (18), shows the dependence of 
J^(l) /j-fo(o) on m - w for a typical value of k. This figure is of particular 
importance because it also relates the critical current of the system (through 
the relative peaking) to the spacing between coplanar films. The general fea- 
tures of the curve are true for all values of d/w considered in this report. 

If m - w > 20d, J b (l)/J b (0) and hence I c remain approximately constant. 

From equations (20) and (21) it is apparent that as m - w is decreased 

/ h \ 

P ' can be greatly reduced, particularly if k ^ k Q , the critical value for 
convergence of equation (16). Calculations show that if |3d = 1 and k = 2.00, 
l("k) decreases fty 15 percent as m - w is decreased from 20d to d. If 

pd = •= and k = 8.00, I' ’ decreases fty 22 percent in this same interval. No 

Zj C 

dependence of l( a ) on q was found, nor should it he expected, since the maxi- 
m-urn peaking occurs on the outside edges of the system and consequently is in- 
sensitive to q. The variation of with m - w can he reduced hy select- 

ing a smaller k (for fixed 0d) or a smaller (3d (for fixed k). 
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1.45i— 


J b (D 

J b (0) 


1.40 


L 35 


1.30 


0 Id 


5d 


lOd 
m - w 


J 

20d 


Figure 7. - Relative peaking of current density at inside edges of films 
in array b plotted against separation distance between adjacent 
films in units of film thickness for $d = 1, k = 1. 5, and d/w < 0. 001. 


Two approximations were made in the development of the solution and require 
brief explanations. The first assumption was that J does not vary signifi- 
cantly with x. It is easily proven that in the central portion of the films 
J ~ cosh p(x - d/2), thus the assumption is not strictly correct. However, it 
is shown in reference 4 that if this variation is substituted into the original 
two-dimensional integral equation, the subsequent one -dimensional integral equa- 
tion is of the same form as equation (9a), with A replaced by aA. The param- 
eter a is bounded, and depends on the thickness of the film in the following 
manner : 


1 


< a < 


sinh 


Pd 

2 


Pd 

2 


( 22 ) 


Since pd < 1, it follows that 1 < a < 1.04. This would result, at most then, 
in a 2 percent decrease for j(u) in the central region of the films. This 
error is less near the edges of the films since in these regions J(x) more 
nearly approaches the constant value Implicit in the original assumption. 

The second assumption was that equation (14) is valid in the edge regions 
of the films. Although figure 4(b) suggests that the approximation would not be 
accurate here, at the very least equation (17) would provide a reasonable trial 
value for any numerical solution of the integral equation. Starting from equa- 
tion (17), numerical solutions obtained with the aid of an IBM 7094II computer 
revealed that the greatest error in the closed-form solution occurs when k ^ k Q . 
In this region it was found that the use of equation (17) can result in an under- 
estimate of I c of as much as 10 percent. For k < 0.65 k c this error is less 

than 5 percent. It should be noted that the percentage decrease in i(^), pre- 

c 

dieted by the closed-form solution, as m - w is varied from 20d to d, agrees 
with that found by numerical analysis to better than 5 percent. 


CONCLUDING REMARKS 

It can be seen from equation (20) and figure 6 that the presence of a 
second strip line alters the current distribution in the original system only 
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in what would he the inside edges of the new system. Consider now the more 
general case of a 2 by n matrix (i.e., n equally spaced films over a super- 
conducting ground plane, or n closely spaced strip lines). If the direction 
of the current alternates from film to film in a given row, the critical current 
of any 2 by 2 submatrix should be the same as that found for an isolated 2 by 2 
matrix. In this case the additional superconductors only serve to bring the 
current peaking on the outside edges of the submatrix up to the same value as 
appears on the inside edges. Since the maximum peaking occurred previously on 
the inside edge, the critical current of the submatrix would be unaltered. On 
the other hand, if the current is in the same direction in each of the coplanar 
films, the peaking on the outside edges would be decreased to the level found 
previously on the inside edges, thus increasing the critical current of the sub- 
matrix. With arguments such as these, equations (20) and (21) can be used to 
determine the critical current for the more general 2 by n array. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, December 6, 1965. 
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